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Abstract 

We investigate the new agegraphic dark energy (NADE) model with power-law corrected 
entropy in the framework of Hofava-Lifshitz cosmology. For a non-flat universe containing the 
interacting power-law entropy-corrected NADE (PLECNADE) with dark matter, we obtain 
the differential equation of the evolution of density parameter as well as the deceleration 
parameter. To study parametric behavior, we use an interesting form of state parameter 
as function of redshift lo\(z) = luq + u)\Z. We find that phantom crossing occurs for the 
state parameter for a non-zero coupling parameter, thus supporting interacting dark energy 
model. 

PACS numbers: 95.36.+X, 04.60.Pp 
Key words: Dark energy, Quantum gravity 



*KKarami@uok.ac.ir 
^sheykhiQmail.uk.ac.ir 
* mj amil@camp .mist .edu .pk 
§ rmyrzakulov@csufresno.edu 



1 



1 Introduction 



Observational data of type la supernovae (SNela) collected by Riess et al. [1] in the High-redshift 
Supernova Search Team and by Perlmutter et al. [2] in the Supernova Cosmology Project 
Team independently reported that the present observable universe is undergoing an accelerated 
expansion phase. The exotic source for this cosmic acceleration is generally dubbed "dark 
energy" (DE) which is distinguished from ordinary matter (such as baryons and radiation), in 
the sense that it has negative pressure. This negative pressure leads to the accelerated expansion 
of the universe by counteracting the gravitational force. The astrophysical observations show 
that about 70% of the present energy of the universe is contained in DE. Although the nature 
and cosmological origin of DE is still enigmatic at the present, a great variety of models has 
been proposed to describe the DE (see e.g., the reviews [3, 4]). Two promising candidates are 
the holographic DE (HDE) [5] and the agegraphic DE (ADE) [6] models which are originated 
from some considerations of the features of the quantum theory of gravity. 

It is curious to note that the HDE model has its origin (i.e. definition and derivation) 
depends on the Bekenstein-Hawking (BH) entropy-area relationship Sbh = A/ 4 of black hole 
thermodynamics, where A is the area of the horizon [7]. However, this definition of HDE can be 
modified (or corrected) due to the various correction procedures applied to gravity theories [8]. 
For instance, the corrections to the entropy which appear in dealing with the entanglement of 
quantum fields in and out the horizon [9] generate a power-corrected area term in the entropy 
expression. The power-law corrected entropy has the form [10] 



s = j 



l-K a A 



(1) 



where a is a dimensionless constant whose value is currently under debate and determining its 
unique and precise value requires separate investigation, and 



a(47r) 
(4 - a)rf 



where r c is the crossover scale. The second term in Eq. (1) can be regarded as a power- 
law correction to the entropy-area law, resulting from entanglement i.e. the wave function of 
the field is taken to be a superposition/entanglement of ground and exited states [9]. The 
entanglement entropy of the ground state satisfies the BH entropy-area relationship. Only the 
excited state contributes to the correction, and more excitations produce more deviation from 
the BH entropy-area law [11] (also see [12] for a review on the origin of black hole entropy 
through entanglement). This lends further credence to entanglement as a possible source of 
black hole entropy. The correction term is also more significant for higher excitations [9]. It is 
important to note that the correction term falls off rapidly with increasing A. So for large black 
holes the correction term falls off rapidly and the BH entropy-area law is recovered, whereas for 
the small black holes the correction is significant. 

The ADE model is originated from the uncertainty relation of quantum mechanics together 
with the gravitational effect in general relativity (GR). The ADE model assumes that the ob- 
served DE comes from the spacetime and matter field fluctuations in the universe [6] . Following 
the line of quantum fluctuations of spacetime, Karolyhazy [13] proposed that the distance in 
Minkowski spacetime cannot be known to a better accuracy than St = where e is a 

dimensionless constant of order unity and tp is the reduced Planck time. Based on Karoly- 
hazy relation, Maziashvili proposed that the energy density of metric fluctuations of Minkowski 
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spacetime is given by [14] 

where Mp is the reduced Planck mass Since in the original ADE model the age 

of the universe is chosen as the length measure, instead of the horizon distance, the causality 
problem in the HDE is avoided [6]. The original ADE model had some difficulties. In particular, 
it cannot justify the matter-dominated era [6]. This motivated Wei and Cai [15] to propose the 
new ADE (NADE) model, while the time scale is chosen to be the conformal time instead of 
the age of the universe. The NADE density is given by [15] 



PA = (4) 



V 

where 3n 2 is the numerical factor and rj is the conformal time and defined as 

fdt f a da 

r] = J^ = J Ha^- (5) 

The ADE models have been examined and constrained by various astronomical observations 
[16, 17, 18, 19, 20, 21]. Inspired by the power-law corrected entropy relation (1), and following 
the derivation of HDE [22] and entropy-corrected HDE (ECHDE) [23], we can easily obtain the 
so-called "power-law entropy-corrected" NADE (PLECNADE) whose the scale is chosen to be 
the conformal time rj. Therefore, we write down the energy density of PLECNADE as [24, 25] 

3n 2 M 2 I3M 2 
PA = — , (6) 

where (3 is a dimensional constant whose the precise value needs to be determined. In this paper, 
our aim is to investigate the PLECNADE model in Hofava-Lifshitz cosmology. 

The plane of the paper as follows: In section 2, we give a brief review of Hofava-Lifshitz 
cosmology in detailed balance case. In section 3, we construct a model of interaction between 
DE and DM. In section 4, we discuss some cosmological implications of this model. We obtain 
the evolution of dimensionless energy density, deceleration parameter and equation of state 
parameter of PLECNADE model. In section 5 we give the conclusion. 



2 Basics of Hofava-Lifshitz Cosmology: Detailed Balance Case 

Recently a power-counting renormalizable UV complete theory of gravity was proposed by 
Hof ava [26] . Quantum gravity models based on an anisotropic scaling of the space and time di- 
mensions have recently attracted significant attention [27]. In particular, Hofava-Lifshitz point 
gravity might be has desirable features, but in its original incarnation one is forced to accept a 
non-zero cosmological constant of the wrong sign to be compatible with observations [28]. At a 
first look it seems that this non-relativistic model for quantum gravity has a well defined IR limit 
and it reduces to GR. But as it was first indicated by Mukohyama [29], Hofava-Lifshitz theory 
mimics GR plus dark matter (DM). This theory has a scale invariant power spectrum which de- 
scribes inflation [30]. Moreover, some new integrable and nonintegrable cosmological models of 
the Hofava-Lifshitz gravity have been discussed in [31]. Phenomenologically, in Hofava-Lifshitz 
gravity the radiation energy density decreases proportional to a~ 6 [32]. Hence the resultant 
baryon asymmetry as well as the stochastic gravity waves can be enhanced. Some cosmological 
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solutions to in Hof ava-Lifshitz gravity are obtained previously [33] . Saridakis formulated Horava- 
Lifshitz cosmology with an additional scalar field and showed that Horava-Lifshitz dark energy 
naturally presents very interesting behaviors, possessing a varying equation-of-state parameter, 
exhibiting phantom behavior and allowing for a realization of the phantom divide crossing. In 
addition, Horava-Lifshitz dark energy guarantees for a bounce at small scale factors and it may 
trigger the turnaround at large scale factors, leading naturally to cyclic cosmology [34]. The 
large scale evolution and curvature perturbations in HL gravity are explored in [35], while the 
origin of primordial large-scale magnetic fields in the Horava's non-relativistic gravity have been 
discussed in [36]. The generalized second law of thermodynamics in Horava-Lifshitz cosmology 
is studied in [37]. For reviews on the scenario where the cosmological evolution is governed by 
Horava-Lifshitz gravity see [29, 38]. 

Under the detailed balance and the project ability conditions, the modified Friedmann equa- 
tions in the framework of Horava-Lifshitz (HL) gravity are given by [26] 



H 



6(3A - 1 



-Pm + 



6(3A - 1) 



3n 2 p 2 k 2 3kVA 2 



+ 



8(3A - l)a 4 8(3A - 1 



K 4 fi 2 Ak 



8(3A- l) 2 a 2 ' 



(7) 



• 3 o K K 

H + -H 2 = -rp 



2 4(3A-ir m 4(3A-1) 



K 2 n 2 k 2 3k 2 h 2 A 2 



8(3A - l)a 4 8(3A - 1) 



K 4 n 2 Ak ^ 



16(3A-l) 2 a 2 ' 



where H = ^ is the Hubble parameter, A is a dimensionless constant and A is a positive constant 
which as usual is related to the cosmological constant in the IR limit. The parameters k and 
fi are constants. Also k denotes the curvature of space k = 0,1, — 1 for a flat, closed an open 
universe, respectively. Furthermore, p m and p m are the energy density and pressure of the 
matter. 

Noticing the form of the above Friedmann equations, we can define the energy density 
and pressure pa for DE as 

3K 2 p 2 k 2 3^VA 2 
PA ~ 8(3A - l)a 4 + 8(3A - 1) ' ( ) 



pa - ,:s:, 4 - ^f^r- (io) 



*Vfc 2 3«V 2 A 2 
8(3A-l)a 4 8(3A-1)' 

The first term on the right hand side proportional to a~ 4 is effectively the "dark radiation 
term" , present in HL cosmology [39] , while the second term is referred as an explicit cosmological 
constant. 

Finally in order for these expressions to match with the standard Friedmann equations we 
set [39, 40] 

k 2 

° c = 16vr(3A-l)' 

KVA =1, (12) 



8(3A- l) 2 



where G c is the "cosmological" Newton's constant. Note that in gravitational theories with 
the violation of Lorentz invariance (like HL gravity) the "gravitational" Newton's constant G g , 
which is present in the gravitational action, differs from the "cosmological" Newton's constant 
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G c , which is present in the Friedmann equations, unless Lorentz invariance is restored [41]. For 
the sake of completeness we write 



G s = 



32vr' 



(13) 



Note that in the IR limit (A = 1), where Lorentz invariance is restored, G c and G g are the same. 
Further we can rewrite the modified Friedmann Eqs. (7) and (8) in the usual forms as 



H 2 + 



3 2 k 



8vrG r 



(pm + PA), 



(14) 
(15) 



3 Model with Interaction 



Here we would like to investigate the PLECNADE in HL theory. To do this we consider a 
spatially non-flat Friedmann-Robertson- Walker (FRW) universe containing the PLECNADE 
and DM. Let us define the dimensionless energy densities as 



a. 



8vrG r 



p A 8vrG c 
Per 3H Z 



Per 3F 2 ™ 

thus the Friedmann Eq. (7) can be rewritten as 

1 — £l k = A + r^- 
Taking time derivative of Eq. (6) and using relation fj = l/a, we get 



fit — 



k 



f>\ ( — 

,ar) 



-2/9A + 



ml 



(a -2) 



(16) 



(17) 



(18) 



Also, if we take the time derivative of the second relation in Eq. (16) after using (18), as well 
as relations rj = l/a and &a = H£l' A , we obtain the equation of motion for U\ as 



l fT 2 aHr] GgSaH 3 ^ 1 



(19) 



Here, prime denotes the derivative with respect to x = In a. Taking derivative of 
—k/(a 2 H 2 ) with respect to x = In a, one gets 



(20) 



To be more general we consider an interaction between DM and PLECNADE. The recent obser- 
vational evidence provided by the galaxy clusters supports the interaction between DE and DM 
[42]. In this case, the energy densities of DE and DM no longer satisfy independent conservation 
laws. They obey instead 

PA + SH(1 + u>\)pa = —Q, (21) 



Pm + 3Fp m = Q, 



(22) 
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where Q = 3b 2 Hp A stands for the interaction term with coupling constant b 2 . Note that the form 
of Q is chosen purely phenomenologically i.e. to obtain certain desirable cosmological findings 
including phantom crossing and accelerated expansion. In literature, one can find numerous 
forms of Q(Hp) while we chose a simpler form sufficient for our purpose. A more general form 
of Q was proposed in [43]. Also the three interacting fluids has been studied before to investigate 
the triple coincidence problem [44]. Differentiating the Friedmann Eq. (7) with respect to time 
and using Eqs. (16), (17), (18), (21) and (22) we find 



H 

H 2 



; \n k - 3(1 - n A ) + 3b 2 n A 

Inserting this result into Eq. (19) one gets 

n' A = n A [3(1 - n A ) - 3& 2 ^a - n fc ] + (1 - n A ) 

Combining Eq. (23) with (20) we have 



n A G c f3(a - 2) 



aHr] G g 6atf V +1 ' 



-2VL A G c P(a - 2) 



aHr] G g 3aH 3 rf a+1 



Q' k = Q k 



(1 - n k ) - zn A - 3b 2 n A + 



2tt A G c f3(a-2) 



aHr] G g 3aif V +1 



Adding Eqs. (24) and (25) yields 



n' A + n' k = (1 - n k - n A ) 



£l k + 3il A 



3b 2 n A (n k + n A ) g c /%-2) 2n A 



+ 



(l-n k -n A ) G g 3a#V +1 aHr] 
For completeness we give the deceleration parameter which is defined as 



After combining Eq. (23) with (27) we get 

q = \[l ~ V k - 3(1 + b 2 )Q A ] + 
Using definitions (16) as well as Eq. (17) we have 



n A G c p(a-2) 



aHr] G g 3a#V +1 ' 



which from it we can obtain 



_ Pm _ Pm 

(1 — fifc — £l A ) 



dln/9 A p' m tt' m tt'j 



Q A , 



d In a 



<>, ' 



£l A 



4 Cosmological Implications 



(23) 



(24) 



(25) 



(26) 



(27) 



(28) 

(29) 
(30) 



In this section, we study some cosmological consequences of a phenomenologically time-dependent 
parameterization for the PLECNADE equation of state as 



ui A (z) = OJq + UJiZ. 



(31) 
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It was shown in [45] that this parameterization allows to divide the parametric plane (coo,ui) in 
defined regions associated to distinct classes of DE models that can be confirmed or excluded 
from a confrontation with current observational data. 

After using Eq. (21), the evolution of the DE density is obtained as [3, 46] 



PA 
PAo 



a -3(l+o) -wi +b 2 ) e 3w! z 



The Taylor expansion of the DE density around do = 1 at the present time yields 

din p A 



In p A = In p Ao + 



d In a 



Id 2 In p A 
In a + -— -7j 

o 2d(lna) 2 



(In a) 2 + 



(32) 



(33) 



Using the fact that for small redshifts, In a = — ln(l + z) 
respectively, reduce to 

Mpa/pAo) on , , tf\ 3 
= -3(1 + uj + 6 Z ) - -uiz, 

In a 2 



ln(p A /p Ao ) _ din pa Id 2 In p A 



o 2d(lna) 2 



-z + %, Eqs. (32) and (33), 

(34) 
(35) 



z. 



In a d In a 

Comparing Eq. (34) with (35), we find that these two equations are consistent provided we have 

1 din pa 



w = - 



3 din a 



-l-b z 



Id 2 In p A 

"l = a 



3d(lna) 

Inserting Eq. (30) in (36) and (37), after using (22), yields 



1 



ui = -- 



n A + n> k 



n A (i-n k - n A ) 



1 — Q fc — J7a / o ' 



(36) 
(37) 

(38) 



"1 = 3 



3b 2 n' A 3b 2 n A (n' A + 
(i-n k - n A ) + (i - n fc - o A ) 2 + sh. ~ 

(iy A + ^,) 2 



"a + «Z 



+ 



(i - n k - n A ) (i - n fc - oo 2 

Substituting Eqs. (24) and (26) into (38) we reach 



J o 



w 



G c (3(a 



2n Ao \ 



1. 



(39) 



3n Ao v g § 3 ^o +1 ^o%; 

Taking derivative of Eqs. (24) and (25) with respect to x = In a and using (40), one gets 

n A = -(n A n k + n f A n k )-3n A (n A -n' A -i)(cj + b 2 + i) + (i-n A )A 

-6Q' A n A (b 2 + 1) + 3Q' A , 

and 



(40) 



(41) 



n'l = n' k (i - 2n k ) - 3n A (n k - n' k )(oj + b 2 + i)- Ati k - 3(n' k n A + n k n' A )(b 2 + 1), (42) 
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where A is given by 



G c P(a - 2) (3H + (a + 1)?? 



aHr]\H 2 Hr] J aHr] G g 3aH 3 r] a+1 \H 2 
The above expression for A can also be rewritten as 



Hr] 



(43) 



A = 



aHr] 

G c (3(a-2 



n k n A - 3fi A (l - «A - & 2 ^a)1 + 3O A (w + 6 2 + 1) 
2(1 + 9) 



aHr] 



2fl A -q-1 



a 



aHr] 



(44) 



where we have used Eqs. (23), (24) as well as the relation r) = 1/a. Adding Eqs. (41) and (42) 
gives 



(i - n A - ft*) [ni + 3fi A (i + w + 6 2 ) + a 

+3f] A [l-(^ fc + ^ A )(l + 6 2 



+ (3f2 A wo-n fc )(n' fc + J2^) 



Finally, by combining Eqs. (24), (26), (41), (44) and (45) with (39) we find 

A 



^ = {1 + ^ + b 2 ) [3io (n A o - 1) - V ko - 3b 2 + 1 
which more explicitly can be written as 



+ 



3fJ Ao ' 



(45) 



(46) 



u>i = (1 + w + b 2 ) 0, Ao (3w - 2) - 3(w + 6 2 ) - O, ko ~ Qo 



+ 



1 



#0% 



/3(a 



G g 9ff 3 ^ +1 fiAo 



a 



- 2(1 + qo) 



(47) 



Therefore, with wo and w\ at hand we can easily write down the explicit expression for oj a (z) in 
Eq. (31) in terms of model parameters such as fi A , Q^, the running parameter A of HL gravity, 
the parameter n of PLECNADE, the interaction coupling b 2 , and the correction coefficients a 
and j3. From Fig. 1 we notice that in the absence of interaction b 2 = 0, ujq = —0.99 showing 
quintessence state. By introducing interaction term, the state parameter evolves to a phantom 
state and gradually goes to more super-phantom state. Figure 1 also shows that the phantom 
crossing for loq happens for b 2 = 0.01 which is compatible with the observation [47]. This is also 
in agreement with the result obtained by [48]. Note that phantom crossing has sound empirical 
support: analysis of Gold SNe and other observational datasets suggests that oj(z) < —1, for 
< z < 0.5 [52]. From Fig. 2, we notice that in the absence of interaction, the first order 
correction u\ to state parameter behaves like quintessence but when interaction is introduced, 
the parameter u\ evolves towards —1, cosmological constant. Here it is probable that oj\ can 
cross the cosmological constant boundary if b 2 > 0.20. For reader's clarity, we did not use 
any initial conditions for plotting both figures since Eq. (40) and Eq. (47) are not differential 
equations. However the behavior of curves in figures is strongly sensitive to the values of free 
parameters. 
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5 Conclusions 



It has been shown that the origin of black hole entropy may lie in the entanglement of quantum 
fields between inside and outside of the horizon [9]. Since the modes of gravitational fluctuations 
in a black hole background behave as scalar fields, one is able to compute the entanglement 
entropy of such a field, by tracing over its degrees of freedom inside a sphere. In this way the 
authors of [9] showed that the black hole entropy is proportional to the area of the sphere when 
the field is in its ground state, but a correction term proportional to a fractional power of area 
results when the field is in a superposition of ground and excited states. For large horizon areas, 
these corrections are relatively small and the BH entropy-area law is recovered. 

Here, we investigated the PLECNADE scenario in the framework of HL gravity. We con- 
sidered an arbitrary spatial local curvature for the background geometry and allowed for an 
interaction between the PLECNADE and DM. We obtained the deceleration parameter as well 
as the differential equation which determines the evolution of the PLECNADE density parame- 
ter. Using a low redshift expansion of the EoS parameter of PLECNADE as oj\(z) = ojq + ojiz, 
we calculated ujq and uj\ as functions of the PLECNADE and curvature density parameters, £l\ 
and Qk respectively, of the running parameter A of HL gravity, of the parameter n of PLEC- 
NADE, of the interaction coupling b 2 , and of the coefficients of correction terms a and (3. It is 
quite interesting to note that phantom crossing for ujq happens for b 2 = 0.01 i.e. a small but non- 
zero interaction parameter. In literature, there are two well-studied ways for phantom crossing: 
via modified gravities including scalar-tensor or Gauss-Bonnet braneworld models [53] or by 
introducing an interaction between fluid dark energy and dark matter. In the later models, its a 
generic feature of interacting dark energy models to have phantom crossing, for instance, using 
different forms of dark energy including Chaplygin gas [54], quintessence [55], new-agegraphic 
and holographic dark energy [20]. In [56], the present authors studied the PLECNADE inter- 
acting with matter in Brans-Dicke gravity and obtained the state parameter behaving to cross 
the phantom divide for small values of coupling parameter b. 
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PLECNADE model in HL gravity 
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Figure 1: The EoS parameter uj of the PLECNADE in HL gravity, Eq. (40), versus the 
interacting coupling parameter b 2 . Auxiliary parameters are: n = 2.716 [16], a = —7.5, (3 = 
-14.8, % = 11 [49], ^A = 0.728, n ko = -0.013 [47], A = 1.02 [50], G c /G g = 2/(3A - 1) = 0.97, 
Hq = 74.2 Km S^ 1 Mpc" 1 [51] and M p 2 = 8irG g = 1. 
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Figure 2: The EoS parameter oj\ of the PLECNADE in HL gravity, Eq. (47), versus the 
interacting coupling parameter b 2 . Auxiliary parameters as in Fig. 1. 



13 



